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Abstract. 

Let N > 1 be a fixed integer and (Ci, . . . , Cjv, Q) a random element of M(d X d, R) N X R d . We 
consider solutions of multivariate smoothing transforms, i.e. random variables R satisfying 

N 



where = denotes equality in distribution, and R,Ri,...,Rn are independent identically dis- 
tributed Revalued random variables, and independent of (Ci, . . . , CjV; Q)- We briefly review 
conditions for the existence of solutions, and then study their asymptotic behaviour. We show 
that under natural conditions, these solutions exhibit heavy tails. Our results also cover the case 
of complex valued weights (Ci, . . . , Cjv). 

1. Introduction 

Let N > 1 be a fixed integer and (Ci, . . . , Cn, Q) a random element of M{d x d, R) N x R d . This 
induces a mapping S : 9H(R d ) — > 9H(R d ) on the space of probability measures on M. d , defined by 



(1.1) 



Six := C C * X < + Q^j . 



where X\, . . . ,Xn are independent identically distributed (iid) with distribution fi, independent of 
the weights and the immigration term (C\, . . . , Cat, Q), and C denotes the law of 2<=i Ci^i + Q- 

On a suitably chosen complete metric subspaces of dJl(R d ), S possesses a unique fixed point. In 
terms of random variables, this means 

N 

(1.2) R = Y,CiRi + Q, 

i=l 

with R, Ri, . . . , i?jv iid and independent of (Ci, . . . , Cat, Q). Uniqueness of R then means of course 
uniqueness of its distribution. We will slightly abuse notation, and call the random variable R above 
a solution of the fixed point equation. In fact, its distribution is a fixed point of S. 

We will also consider a particular case of (1.2) when Q = 0, i.e. we will study solutions of the 
equation 

N 



(1.3) R = Y,CiRi 



i=l 
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It turns out that there are some subtle differences between those two cases. To distinguish between 
them we will call the stochastic equation (1.2) inhomogeneous and (1.3) homogeneous. 

The case of nonnegative scalar weights, i.e. Cj £ R+, known as the smoothing transform, has 
drawn much attention, see e.g. the classical works [20, 27, 28], as well as [3] and the references 
therein. Also the case of real-valued scalar weights has been studied, see [17] or very recently in 
[29]. 

The study of its multivariate analogue has a much shorter history, but it draws more and more 
attention. We will focus on two different kinds of assumptions, the first being Ci from the set of 
similarities, i.e. products of orthogonal matrices and dilations. This particularly covers the case of 
the smoothing transform with complex valued weights. An equation of this type was very recently 
studied in [18]. Another situation, where similarities appear as weights, is the joint distribution of 
two statistics appearing in phylogenetic trees, see [10, Equation (14)]; or the joint distribution of 
key comparisons and key exchanges for Quicksort, see [33, Theorem 4.1]. In all these papers, only 
the case of solutions with exponential moments has been studied. In particular in the light of the 
article [25], it is tempting to search for solutions with finite expectation, but heavy tails. 

Secondly, we will treat general matrices, under some density assumptions on their distribution, 
as introduced in [4]. There, new properties, unknown from the one-dimensional case as well as prior 
multidimensional studies concerning positive solutions ([13, 30]), appear. We will discuss them 
by a nice example, which also allows us to (partially) answer a question raised by Neininger and 
Ruschendorff [33, Problem 3.2]. 

A particular motivation for our work are equilibrium distributions of kinetic models. They can 
be studied as fixed points of smoothing transforms, as discussed for the one-dimensional case in the 
very recent article [8]. It is stated there [8, Theorem 2.2 (iii, iv)] that these distributions show heavy 
tail behaviour in the sense that moments of high order are infinite. Our results (for similarities) 
then give precise tail asymptotics for these models. We will describe the basic ideas in subsection 
8.1, and derive heuristically a stochastic fixed point equation for the distribution of particle velocity 
in Maxwell gas. Naturally, this is a distribution on M. d , and it was studied in [9, Example 6.1], where 
it was shown to have heavy tails, [9, 6.1.4], but without precise tail index. Our results in the case 
of general matrices give these tail asympotics. The details are worked out in subsection 8.1. 

1.1. What can be expected from the one-dimensional case? The main properties of the 
fixed points are governed by the function rh(s) := E^ili > which is convex with m(0) = N > 1, 
thus there are at most two solutions < a < f3 to rh(s) = 1. It is shown in [3, Theorem 8.1], that 
the solutions to (1.2) are of the form 

X = R + hW 1/a Y, 

where ft > and Y is a one-sided stable law of index a independent of (R, W), if a € (0, 1), and 
Y = 1 if a = 1. The random variables (i?, W) can be expressed in terms of the weighted branching 
process (WBP) associated with (Ci, . . . ,Cn,Q)- The WBP will be defined in Section 3. Note 
that in the recent preprint [29], there is a similar result for real- valued weights, and vector valued 
solutions. This corresponds to diagonal matrices (which are of course similarities) in our setting. 

If ft, > 0, tails of the solutions are governed by a. If ft, = 0, then the solution is given by R, and 
it is this special solution (called minimal solution in [3]), which we are interested in. It was first 
shown by Guivarc'h [22] (for the homogeneous case (1.3) and a = 1) that given the existence of /3, 
the solution R has heavy tails with index j3, i.e. 

(1.4) lim t^P (R>t)= K. 

t— ¥oa 

For the inhomogeneous equation, and also for the case of real-valued weights, this was recently 
shown by Jelenkovic and Olvera-Cravioto [26, 25] (see also [15, 16]). 
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We extend their work to the multivariate case, and give in particular a nice holomorphic argument 
that shows K > 0, at least for similarities, in the general case it gives an equivalent condition. The 
main part of this work is devoted to the study of the case of general regular matrices; most of the 
proofs carry over to the case of similarities, which will be treated separately. 

1.2. Acknowledgements. The authors want to thank Gerold Alsmeyer and Jacek Zienkiewicz for 
helpful discussions during the preparation of the paper, and Daniel Matthes for information about 
kinetic models and pointing out reference [9]. We are grateful to the referee for a very careful reading 
of the manuscript and many helpful suggestions improving the presentation. 



2. Statement of results 

2.1. Notation and assumptions. The Euclidean space Mr is endowed with the scalar product 
(x,y) = ^2i=i x iUi ana - the norm \x\ = \J (x, x). The unit sphere in R d is denoted by S, and the 
projection of a vector x <E K d \ {0} on S by (x)~ := \x\ 1 x. Moreover, on the space M(d x d, M) of 
dx d- real matrices we will consider the operator norm || • || associated with the Euclidean norm | • 
on M. d , i.e. \\a\\ = sup xeS \ax\ for every dx d matrix a. The open ball of radius 8 around a matrix A 
w.r.t. to the operator norm is denoted by B$(A). The Lebesgue measure on M(d x d,M), identified 
with R d \ is denoted by X dxd . We abbreviate R+ = (0, oo). 

The set of continuous functions / : E — > R on a metric space E is denoted by C(E), and by 
Cc{E) we denote the compactly supported continuous functions / : E — >• M. The set of fc-times 
(Frechet) differentiable functions is denoted by C k (E). 

Let J be a uniformly distributed random variable on the set {1,...,A}, independent of the 

random variable (Ci, . . . , Cm, Q)- We introduce the random variable C = Ci, together with a 
sequence (CW) n > of iid copies of C. Products of . . . , CW will be denoted by LT„ := C (1) ■ 
. . . • . We will see that the right multivariate expression for to is given by the function 

(2.1) m(s) :=N lim (E||n„]| s )" , 

which is defined for < s < Soo := sup{s > : E ||C|| S < oo}. Note that Soo = sup{s > : 
maxE ||Ci|| s < oo}. Again, to is a convex function, and we define 

(2.2) a := inf{s > : m(s) < 1} 

(2.3) 13 := sup{s > : m(s) < 1}. 

We will always assume that a < j3 < Sqc, then m(a) = m((3) = 1. Though interesting examples with 
/3 = Sqo exist, in applications, to usually has to be approximated by simulations, so it is convenient 
to introduce the standing assumption /3 < Soo for it will also be needed in parts of the proofs. Later 
on, we will see that to is differentiable (at /3), then m'(/3) > 0. 

Applying the random variable Cj one can easily prove (Proposition A.l) that Ci, . . . , Cjy are 
identically distributed (but dependent), which we w.l.o.g. assume from now on. 

It may happen that the distribution of R is degenerate. We exclude this case by the assumption 

(not triviality) For all r E R d , P b- = ^ C k r + Q^j < 1. 

As soon as the first moments of R, (Ci, . . . , Cn,Q) exist (1.2) yields an identity for the expectation 
of R. In some cases, we have to assume that this identity indeed has a solution, i.e. there is r £ R d 
such that 



(eigenvalue) 



r = ATECY + EQ. 



4 



D. BURACZEWSKI ET AL. 



Moreover, if second moments exists, the covariance matrix S of a solution R to the homogeneous 
equation (1.3) has to satisfy 

(covariance) £ = N E[C*£C T ] . 

2.2. Existence and uniqueness results. The following existence and uniqueness results are ob- 
tained from results of [32]. For the reader's convenience we give some ideas in Appendix B. 

Proposition 2.4 (homogeneous case). Let (Ci,...,Cjv) be a random element of M(d x d, M.) N . 
Assume a < min{2,Soo} and let r £ R d , satisfy (eigenvalue). 

• Case a < 2 : Then there is a unique solution R to the homogeneous equation (1.3) with 
Ei? = r, and E \R\ S < oo for all s < min{soo,/3}. R = if and only if r = . 

• Case a = 2: Assume in addition that there is a symmetric and positive definite matrix £ £ 
M(d x d, R), satisfying (covariance). Then there is a unique solution R to the homogeneous 
equation (1.3) with MR = r, covariance matrix K(i?) = £ and E|i?| s < oo /or all s < 
min{s 00 ,^}. 

Proposition 2.5 (inhomogeneous case). Let (Ci, . . . , Cat, Q) be a random element of M(d x 
d,R) N x R d . Assume < a < minjsoo^}, and, if a > 1, let r £ R d satisfy (eigenvalue). If 
E|Q| S < oo for some a < s < minjsoo^} then there is a unique solution R to the inhomogeneous 
equation (1.2) with Ei? = r. 

Moreover if E|Q| min { s ~^> < oo, then E\R\ S < oo for all s < min{ Soo ,/3}. 

Remark 2.1. In dimension d = 1, exploiting the identity (covariance) for the variance of a solution 
R, one can show that if a > 2, then there is no nontrivial solution with finite moment of order 
s > a. See [2] for a detailed discussion. Without further assumptions, this is not true in dimension 
d > 2, as the following (even deterministic) example shows: Let 

m (V A) ■ 

Then m(s) = A^Cif = N(N~^) S = N^i , thus a = 3. But if R 2 has standard normal distri- 
bution, the random vector (0, R 2 ) T is obviously a solution of the smoothing transform associated 
with (2.6). The point is that m is only concerned with the largest eigenvalue of C, but there may 
be solutions concentrated on subspaces. 

These existence and uniqueness results hold for general d x d-matrices {C\ , . . . , Cn ) £ M(d X d, R) . 
In order to describe the asymptotic behavior of R, we will also consider more specialised classes of 
matrices, namely invertible matrices, i.e. (Ci, . . . , Cat) £ GL(d, M), and the subclass of similarities. 

The limit theorems are now stated separately for two main cases: matrices satisfying some 
irreducible and density hypotheses (here two subcases appear), and matrices being similarities. 

2.3. Asymptotic behavior - general matrices. We start with the first case studied recently in 
[4]. We will assume that C acts irreducibly on the sphere, i.e. 

(irred) V xeS V ope „ c/cs maxP ((aT^T G U) > 0, 

n>l 

and that the law of C is spread out, i.e. 

(density) ^AeGL(d,R) ^c,s>o 3„ eN P (n„ £ •) > cl Bs{A) \ dxd . 
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Theorem 2.7. Let (Ci, . . . , Cn, Q) be a random element of GL(d, R) x R d , and R be the unique 
solution to (1.2) given by Propositions 2.4, resp. 2.5. Let the assumptions of these propositions hold, 
and assume the existence of a < (3 < such that m(a) = m((3) = 1, as well as E\Q\ Sx ' < oo. 
Let the conditions (irred),(density) hold. 
Then 

(2.8) lim t^F(xR > t) = K- e(x), 

t— >oc 

for a positive continuous function e : S i— > (0, oo). K > if and only i/E \R\ = oo. 

In the light of [7, 11], this property is (as long as f3 is not an even integer) equivalent to multivariate 
regular variation: There is a unique Radon measure A on R d , such that for every / G Cc(R d \ {0}), 
the set of compactly supported functions, 

lim t^Ef(t- 1 R) = [ f(x)A(dx). 

JR d \{0} 

Particularly motivated by [9, Example 6.1], we additionally state the theorem with a slightly 
different set of assumptions, which are in particular suitable for random matrices of the type C = 
YY T , where Y is a random vector with a spread-out density. 

So here we do not assume C £ GL(d,M), instead we need: 

(not vanishing) \/ x£S P (xC = 0) = 0. 

A probability measure on K is called nonarithmetic, if the smallest closed group containing its 
support is R. The densitiy assumption now reads as follows: 

(minorization) 3 p>0 3 0eOT(s) 3 v , eOT(K)> nonarit hmetic VxeS P {(xO)~ € •, log \xC\ G ■) > p <j> ® ip. 

Theorem 2.9. Let {C\, . . . , Cn, Q) be a random element of M(d xd,R) N xR d , and R be the unique 
solution to (1.2) given by Propositions 2.4, resp. 2.5. Let the assumptions of these propositions hold, 
and assume the existence of a < (3 < Soo such that m{a) = m(/3) = 1, as well as E |Q| Soc < oo. Let 
the conditions (not vanishing),(minorization) hold. 
Then 

(2.10) lim t F(xR>t) = K-e(x), 

t— >oo 

for a positive continuous function e : S i— > (0, oo). Moreover, K > if and only j/E \R\ = oo. 

2.4. Asymptotic behavior - group of similarities. Next we will assume that Ci are elements 
of the group of similarities i.e. the group of elements g of GL(d, R) satisfying 

\gx\ = \\g\\\x\ 

for every xelf The group of similarities G is the direct product of the multiplicative group R + 
and the orthogonal group O(d). 

Notice that if g G G then its norm is given by its radial part, i.e. if g = tk, t G R + , k G 0(d), 
then ||<7|| = t. This implies in particular that in this case m(s) = ATE||C|| S . We denote by /x the law 
of C and by G p the subgroup of G generated by the support of /i. We will assume that log ||C|| is 
nonarithmetic. Then there is a subgroup of 0(d) such that 

Gfi = R + x Hp. 

(see [14] for more about the structure of the group G). 
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Theorem 2.11. Let (Ci, . . . , Cjv, Q) be a random element ofG N xM. d , and R be the unique solution 
to (1.2) given by Propositions 2.4, resp. 2.5. Let the assumptions of these propositions hold, and 
assume the existence of a < (3 < Soo such thatm(a) — m(j3) = 1. Assume moreover thatK\Q\^ < oo 
and that log || C'i | are nonarithmetic. Then 

lim ^P(|U| >t)=K + . 



t— > oo 

P d 



Moreover there exists a unique Radon measure A such that for any f £ Cc(IR \ {0}) 

lim \\a\\~ P Ef(aR)= [ f(x)A(dx). 



||a||->0, aeG„ JK d \{0} 

The measure A is homogeneous and A(dg) = cr(dk)j0^ for some finite measure on S such that 

N /3 N 

CiRi + Q -J2\CiRi^ 



ct(5) = — E 



i=l 



i=l 



where mp = E[\\C\f log \\C\\] > 0. 

Remark 2.1. The theorem stated above is an analogy with the one obtained in [14] for the solution 
of R = CR + Q. 

In this case of similarities, we obtain a much stronger dichotomy concerning nontriviality of the 
limit measure: 

Proposition 2.12. Suppose that the assumptions of Theorem 2.11 are satisfied and additionally 
E|Q| Soc < oo. Then either K + and o~{S) are positive, or E |i?| s < oo for all s < s^. 

From this we will deduce a sufficient condition for er(S') > 0: 

Proposition 2.13. Under the assumptions of Theorem 2.11, let 7 6 (f3, Sqo), be such that E ||C|| 7 = 
1, andE\Q\ s °° < 00. Then 

(2.14) a(S) > <^> (not triviality) holds. 

2.5. The structure of the paper. The organization of the paper is as follows. We start with a 
brief section introducing the weighted branching process, a stochastic model that allows to study 
iterations of S in terms of random variables. The proof of existence and uniqueness of solutions 
to (1.2) and (1.3) (Proposition 2.4 and Proposition 2.5) is classical and we present only main 
arguments, therefore we postpone the proof to Appendix B. We consider separately general matrices 
and similarities. Although ideas are similar, the two proofs require slightly different techniques. All 
proofs are divided into two parts. First we prove the existence of the limit and then the positivity 
of the limiting constant. General matrices are treated in Sections 4 and 5. First we introduce the 
tools (transfer operators) in Section 4 and then in Section 5 we formulate the renewal theorem and 
conclude existence of the limit. Asymptotic behavior of solutions in the similarity case is considered 
in Section 6. Finally we finish proofs of our main results in Section 7, where we prove positivity 
of the limiting constants. Since the proofs consist of tedious calculations, for a better stream of 
arguments, some proofs have been carried over into Appendix C and Appendix D. In Section 8 we 
present applications of our results to very concrete models. 

3. The weighted branching process 

In this section, we introduce a stochastic model, the weighted branching process (WBP), that 
produces a sequence of random variables (Y n ) with C(Y n ) = S n (fj,). We consider random variables, 
indexed by a A^-ary tree T. For a node v = (ii, . . . , if.), we denote its level by = k, its 
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ancestor in the l-th level, I < k, by v\l = (ii, . . . , ii) and its ith successor by vi = (ii, . . . , i). 
The root is denoted by 0. To define the weighted branching process, we assign to each node 
v a copy (Ci(i>), . . . , Cjv(i>), Q) of the vector (Ci, . . . , Cjv, Q), independent of all other random 
variables. The product along a path is defined recursively by L(0) := Id, the identity matrix, and 
L(vi) := L(v)Ci(v). Given a random variable X, we assign to each node also a copy X(v) of X, 
independent of all other random variables. Then 

Tl-l 

Y n := L(v)X(v) + J2 E L ( v ^ 

\v\=n k=Q \v\=k 

is called the weighted branching process associated with X, (Ci, . . . , CV, Q). It is easy to see that 
if X has distribution /i, then S n (fi) = C(Y n ). Moreover, the weights (L(v))i v i— n are dependent, but 
identically distributed with the same distribution as II„. 

In many cases (see e.g. [13, 25, 26, 30]) one can prove that the sequence Y n converges pointwise 
and its limit provides a solution of the smoothing transform. However in order to prove existence 
theorems in full generality, one has to consider S as an operator on an appropriate complete metric 
space and apply the Banach fixed point theorem. The details are worked out in Appendix B. 

4. Transfer operators and change of measure 

This section and the next section are devoted to the partial proof of our main results, Theorems 
2.7 and 2.9. We prove existence of the limit, and the positivity of the constant is postponed to 
Section 7. Thus we assume that the CVs are in GL(d,M) and satisfy (irred), (density) or that the 
Cj's are in M(d x d,M) and satisfy (minorization) and (not vanishing). 

First, we study a family of transfer operators introduced by the action of C on the sphere S, 
and derive harmonic functions, which allow for a change of measure. The results of this section are 
mainly technical, however this is the main tool that will be used in the next section when applying 
the Markov Renewal Theorem. As a by-product, we prove holomorphicity of m(s). 

4.1. Transfer operators. For < s < Sqo we introduce the family of operators on continuous 
functions on the sphere, T s : C(S) —> C(S), defined by 

T s f(x) :=E[f((xC)~)\xC\']. 

They are well defined as mappings on C(S) since C € GL(d,M) resp. (not vanishing) holds. One 
can easily check that powers of these operators are given by 

If f{x) :=E[/((xn„r)|xn„| s ]. 

Observe that T is a Markov transition operator on S, with associated transition kernel defined by 
P(x, A) := P((xIIx)~ £ A) for x <E S and measurable Ac S. For compact subsets D of GL(d,R), we 
further define the substochastic kernels Pd(x, •) := P((xlli)~ 6 IIj £ D). Then from assumptions 
(irred) and (density), we may derive a property very similar to (minorization): 

Lemma 4.1. Suppose (irred) and (density). There is n e N, p > 0, a compact subset D of GL(d,M.) 
and a probability measure with supp(0) = S such that the minorization condition 

(MC) P n (y,-)>P r ](y,-)> P -cf> 

holds for all y G S . 

Proof. A weaker version of the lemma was proven in [4] (Lemma 2.1): there exists for each x G S 
a compact subset D x of GL{d, M) and 5 x ,p x > 0, n x G N, and a probability measure <p x with 
supp(0a;) = Bs x (x) such that 

(4-2) P n *(y,-)>PZl(y,-)>p x <l> x 
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for all y G S . 

Therefore, for each x G S there is an open ball Bg^ and a probability measure <p on it. This 
is an open covering of the compact set S, choose a finite covering around points x±, . . . ,Xk- Set 
n := LCM(n xl , . . . , n Xk ), and define 



L := . G {0, n Xl ,...,n Xk } n : J^Wj^n, w„ 7^ Oj . 

For ui G L, set 

A, := {M x • ... • Af„ : Afj G D Xui ,. . . , M n G D X „ B }, 
with Z)o containing only the identity matrix, and D := lj wei D u . Observe that this union is finite, 
so D is still compact. With p := 1 define p u := n™=i Pz; > an d p := min^gip^ > 0. Finally, 
with := fc _1 (0! + • • • + 0fc), we have 

p%(y, ■) > 4, E ^ ® • • • ® ^:„„ (y, •) > E |7^^ ^ w>- 

□ 

Remark 4.1. It is a direct consequence of the existence and uniqueness results obtained in Appendix 
B that we may as well study iterated versions of the fixed point equation, i.e. 

R= E Hv)R(v) + Yl E L (v)Q(v), 

\v\—n k—0 \v\— k 

where the matrices (L{v))\ v \ =n are not independent, but identically distributed, with the same 
distribution as Yl n , and (R(v)) v ^j- are iid with the same distribution as R and independent of 
(L(v),Q(v)) V £f. Thus, we may w.l.o.g. assume that if (irred) and (density) hold, then 

(MC) holds with n = 1 
- which we will do from now on. 

This result now allows us to deduce the following properties of T s from the results in [4]: 

Proposition 4.3. Let conditions { (irred), (density)^ or { (not vanishing), (minorization)^ hold, and 
< s < Soo. Then the spectral radius and the dominant eigenvalue ofT s are equal to k(s) := jj-m(s). 
There is a unique strictly positive continuous function e s : S 1— » (0, 00) and a unique probability 
measure v s with supp(f s ) = S such that 



T s v s = k(s)u„, T s e s = n(s)e s , / e s (x)v(dx) = I. 

Js 

Moreover, e s (x) = e s (— x) for all x G S. 

Proof. By the minorization condition (MC) with n = 1, there is a compact subset D, and <p,p such 
that 

T s f(x)>E(f((xCr\xC\ s )l D (C)) 

> iam{\xg\ s : x G S, g G D, s G /} • E (f((xC)~ 1 D (C)) 

f(y)P D (x,dy)>cp f f(y)<f>(dy). 
Js 

The same follows from (minorization). This shows that the operator T s is strictly positive, i.e. it 
maps a nonzero nonnegative continuous function to a strictly positive function on S. Now [4, Lemma 
5.2] states the existence of a dominant eigenvalue k(s), equal to the spectral radius; as well as the 
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existence of the corresponding eigenmeasure. In [4, Lemma 5.3], it is shown that k(s) = -^m(s). 
Properties of the eigenfunction e s are shown in [4, Lemma 5.4]. The assertion about the support of 
v s is a direct consequence of the above calculation, since for all / £ C(S), f > 0, / ^ 0, 

K (s) j f{x)v s {dx) = / {T s f){x)v s {dx)>cp f f(y)cf>(dy)>0, 

J S J S J s 

for supp(</>) = S. □ 

As an observation, we note that along the same lines, one can also obtain a uniform minorization 
for the eigenmeasures v a : 

Lemma 4.4. Let I be a com/pact subset of [0,Soo). Then there is p > and a probability measure 
</> with support S, such that v s > p ■ <p for s £ I . In particular, supp(iAs) = S for all s £ [0, Soo)- 

4.2. Markov Random Walks and Change of Measure. A main ingredient in the proofs in 
[25, 26] was application of a suitable renewal theorem. We will use the Markov renewal theorem 
of [1] (see Section 5), which deals with Markov Random Walks: Let {X n ,U n ) n >o be a temporally 
homogeneous Markov chain onSxM such that 

P((X n+1 ,U n+1 ) £ Ax B\X n ,U n ) = P(X n ,Ax B) a.s. 

for all n > and a transition kernel P from S to S x R. Then the associated sequence (X n , V n ) n >o 
with V n = V n -i + U n is also a Markov chain and called the Markov Random Walk with driving 
chain {X n ) n >0- 

We will study the Markov Random Walk defined by 

(x n ,v n ) := ((x n n r,iog|A n„|) = (pr n _iC)~,K-i +io g |x n _ 1 c|). 

Its initial distribution is noted via F x (Xo = x, Vq = 0) = 1. 
Notice that 



(4.5) p Pf(x,t) := — rr-^-Ei (f(X lt Ui + t)e p {X 1 )e^ +t )) 

eiq(x)e pt V / 



is a Markov transition kernel, since 

fpi(x,t) = N E x (l SxR (X 1 ,C/ 1 +t)ep{X 1 )e^ +t A = JL-T p ep{x) = 1. 

The associated probability measure on the path space, and the related expectation E x are 
defined by 

?E X (J(X ,Vo,Xi,V 1 ,...,X n ,V n )) 

(4.6) := —-E x (e^X n )e l3V -f (X ,V ,X 1 ,Vi,...,X n ,V n )), 

ep{x) V J 

for all bounded continuous functions / and all n > 0. 

The transition operator of the driving chain (X n ) n >o under is then given by 

(4-7) epf(x) = ^-(T p e p f)(x). 

From the minorization condition (MC) with n = 1, resp. (minorization) one obtains directly that it 
satisfies the Doeblin condition, i.e. for all x £ S 

P(x,-)>c- 

for the probability measure <f> as defined in (MC). Thus it is in particular positive Harris recurrent. 
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By the regeneration procedure of Athreya and Ney [6], there is a sequence of random times 
(on)n>i> called regeneration epochs, such that for each k > 1, (X CTfc +„) n >o is independent of 
(Xj)o<j<cr k -i with distribution P^ ((X n ) n >o G ■) (see [4, Lemma 4.1]), with <j> defined in Corol- 
lary 4.1 resp. given by (minorization). 

The following lemma is a consequence of the minorization in (MC) via the substochastic kernels 
Pd- 

Lemma 4.8 ([4, Lemma 5.6]). We can choose a sequence of regeneration epochs {o~ n ) n >o, such that 
there is a finite interval I C M with U an G / P-a.s. for all n > 1. 

We then say that V„ has bounded increments at regeneration epochs. 

Remark 4.2. Under (minorization), we even have the stronger property that (X ak+n , f7<7 fc + n )n>o are 
independent of (Xj, Uj)Q<j< ak -i and have distribution P^®^ {{X n , U n ) n >o G ■) . This strong type of 
regeneration is studied in [5]. Obviously, it implies that V n has bounded increments at regeneration 
epochs. 

All in all, we have proven the following result, where the last assertion about tt is a direct 
consequence of the formula (4.7) for 13 P. 

Proposition 4.9. Under the measure ^P as defined in (4.6), (X n ,V n ) n >o is a Markov Random 
Walk, its driving chain (X n ) is a Doeblin chain, thus positive Harris recurrent, and V n has bounded 
increments at regeneration epochs. The stationary distribution -k of (X n ) is given by 

ir(dx) = ep(x)vp(dx). 

4.3. Properties of the operator T z . The final part of this section investigates further the map- 
pings s i — ^ k(s) as well as s i— )• e s and sh>c s . By means of a perturbation theorem for quasicompact 
operators, we will see that these mappings are holomorphic on some small ball around /3. First, we 
show that s H>T S is in fact a holomorphic mapping. For < $tz < Sqq, define the operator T z by 

T z f(x) :=E(\xC\ z f((xCD). 

A family of linear operators T z : B — > B is called weakly holomorphic, if for every x G B, y G B 
{B 1 is the dual space of B) the function 

z h> y{T z x), 

is holomorphic. In our case B = C(S), this is equivalent to being (strongly) holomorphic, see e.g. 
[31, Exercise 8.E.]. 

Lemma 4.10. The mapping z i-> T z is (strongly) holomorphic on the domain < < s^. 

Proof. The proof that T z is holomorphic is taken from [23, p. 253]: Let 7 be a closed path in the 
domain < 5Rz < Sqo, then for any bounded continuous function / and finite measure v, we show 
that f^(v,T z f)dz = 0: 

[ {u, T z f)dz =111 f((Mx)~) \\Mx\\ z P (Cj G dM) v{dx)dz 

Jj Jj JS JGL(d,R) 

f((Mx)~) ( [ \\Mxf dz) P (Cj G dM) v{dx), 

and the innermost integral is zero since z 1— > ||Afa;|| 2 is holomorphic. 

This, together with the fact that z n> (v,T z f) is continuous, implies that T z is weakly holomor- 
phic, thus already strongly holomorphic. □ 
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For a linear operator T denote its spectral radius by r(T). A linear operator T on C(S) is said 
to be quasi- compact, if C(S) can be decomposed into two T-invariant closed subspaces 

C(S) = F®G, 

where r(Tigr) < r(T), while dimF < oo and each eigenvalue of Q\p has modulus r(T). 
Lemma 4.11. The operator Tp is quasi- compact. 

Proof. With the results of Proposition 4.3, what remains to show is the spectral gap property. But 
this follows, since @P satisfies the Doeblin condition, and its spectral properties are one-to-one (resp. 
1-to--^-) with those of Tg, since for all / € C{S), 



T p f{x) = %LPp{f/e p ){x). 



□ 



Then we may apply the perturbation theorem [24, Theorem III. 8] and derive for our situation 
the following corollary: 

Corollary 4.12. There is 5 > 0, such that for all z £ Bg(/3), T z has a simple dominating eigenvalue 
k(z), with eigenfunction e z and eigenmeasure v z , and all mappings 

k:B s (/3)^C, e. : B s [fi) -> C(S) and v. : B s (/3) -> C(S)' 

are holomorphic. In particular, m = Nk is differ entiable in /3. 

5. Implicit Markov Renewal theory on trees 

We will show the following proposition, which gives proof of existence of the limit in Theorems 
2.7 and 2.9: 

Proposition 5.1. Under the assumptions of Theorem 2.7, resp. Theorem 2.9, 



(5.2) lim iP P (xR >t) = f E 

t^oo 2pip Jg 



N 



vCiiU + yQ 



N 



2Ph 

for all x G S, where 

h = (ep((yCD \yC flog \ y c\) Mdy) > o. 

But first, we have to introduce the Markov Renewal Theorem, which we will use in the subsequent 
proof. 

5.1. Markov Renewal Theory. A measurable function g : SxR — > R is called ir-directly Riemann 
integrable if 

(5.3) g{x, ■) is A-a.e. continuous for 7r-almost all x € S 

(5.4) and / sup \g(x,t)\ir(dx) < oo for some i5 > 0, 

Js nez te[nS,(n+l)S) 

where A denotes the Lebesgue measure on R. The following Markov renewal theorem (MRT) is the 
main result of [1]: 
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Theorem 5.5. Let (X n ,V n ) n >o be a nonarithmetic MRW with positive Harris recurrent driving 
chain (X n ) n >o with stationary distribution n. Let I := E^Vi > 0. // g : S x R — > K is a n-directly 
Riemann integrable function, then for ir-almost all i£S, 



(5.6) 



, n>0 



g(u, v) dv n(du). 



S JR 



Remark 5.1. The following extension of the above result is given in [4, Section 7]: If (X n ) n >o is 
a Doeblin chain with bounded increments at regeneration epochs, then the assertion (5.6) holds for 
all x £ S. In our case, this is assured by Proposition 4.9 resp. the subsequent remark. 

For the precise formulation of nonarithmeticity in the context of Markov Random Walks, see 
[35]. We only note here that if (irred), (density) hold, an adaptation of [4, Lemma 5.8] shows that 
(X n , V n )n>o is nonarithmetic under ^P^; and that it is a direct consequence of (minorization). Using 
Proposition (4.9) and the definition of in (4.6), we compute 

'XrVa = J epiyr'E (e ((yCD \yCf log \yC\) ir(dy) 

= J e (ep((yC)~) |i/Cf log|yC|) v p {dy) = l, 3 . 
Using m'(/3) > and [4, Lemma 5.9] yield that lp > 0. 

5.2. Implicit Markov Renewal Theory on Trees. Now we give the proof of Proposition 5.1. 
Let us define 



f(x,t) 



ep(x) 



(xR > e £ 



We would like to write the function / as a potential of some function g and then to apply Theorem 
5.5. However the function / is not sufficiently smooth to satisfy all the hypotheses of the renewal 
theorem. Therefore we consider its smoothed version, i.e. for any function g : S x R — > R we define 
its exponential smoothing 



g{y,t) 



,-(*-») 



g{y,s)ds. 



By [21, Lemma 9.3], if one of f(x,t) and f(x,t), converges for t — > oo, then both of them converge 
to the same limit. So it is sufficient to consider the exponential smoothed version of /. 

Lemma 5.7. The function f satisfies f{x,t) = g*U x (t), where 



g{x,t) 



ep{x) 



(xR > 



NP (xCR > e*) 



Proof. First we expand P (xR > e') into a telescoping sum, using the WBP: 



(xR > 



E 



fc=0 



/ N 

Y p (xL(v)R(v) > e*) - (x 
\v\=k V i=l 



L(yi)R(yi) > e') 



Y P (xL(v)R(v) > e*) 



n-l 



fe=0 



E 

i«i=fc 



N 



vU k R > e*) - Y p (xU k CR > 



N n P (xILnR > e l 
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n-1 

= Nk [ p ( x ke Vk R > e*) - NV (X k e Vk CR > e*)] + N n P (xH n R > e*) 

fe=0 

n— 1 /. 

= J2 Nk v (y R> et ~ v ) - NP (y° R > e *~ v ) p - G ^ e dv ) 

fe=0 ^ 

+iV n P (xn„i? > e*) . 
Multiplying both sides by /ep{x) and using (4.6), we obtain 

/(*, t) = V / — tt" [P (yi? > e*"*) - iVP (j/Ci? > e*- v )] 7V fc ^ , P a (X fc G dy, % G dv) 
-N n V (xIL n R > e*l 



fe=0 



e /Jt 



ep{x)' 

'°Z} f p P(t-v) 

Y / —ty t p > ^ - NF ( yCR > e *~' u )] F * ( Xk e dy > Vk e dv) 



k=0 

+N n —-P (xU n R > e ( ) 

^} r e m 

= Y / 9(V, t - «) ^ (X k G dy, V k G d«) + TV' 1 — -P (a;n„i? > e 4 ) . 
Applying to both sides the exponential smoothing we have 

r ft [3s 

(5.8) f(x, t) = Y 9(V, t - v) "F x (X k G dy, V k G dv) + / e - (t " s) — -7V"P (xiL.i? > e s ) da. 
k=o J J -°° 6/3 W 

Now we want to pass with n to infinity and prove that the second term vanishes. To this purpose 

choose 5 > and uq G N such that N (E ||n n ||^ _ j < 1 — e for all n > no and some e > 0. This 



is possible since m(/3 — <5) < 1. Then, using the Markov inequality, 

e -(*--)_^jVnp( a .n nJ R > e s )ds < / — — iV"P(|n„i?| > e s )ds 

-oo e p( x ) J-oo ep{x) 

ft e -t+(p+i)s E (\Tl n Rf~ S ) ft (s+i)s . 



e/3(x) e 3 ^ 5 ) 7_ 00 ep(x) 

<e st -(l-ey\ 

which tends to zero for each fixed t, as n — > oo. 

The first term in (5.8) for n — > co is almost a renewal function, we want to switch summation 
and integration, to read 

/oo 
g{y,t-v)Y P ^ (Xk G dy, V k G dv) . 

k=0 

Therefore it remains to prove that g is a directly Riemann integrable function. 
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In view of [21] Lemma 9.2, it is sufficient to show sup y6S J R \g(y,s)\X(ds) < oo, i.e. (after a 
change of variables) 



(5.10) 



sup 

yes Jo 



t?- 1 \F (yR >t)-N¥ {yTR >t)\dt<oo. 



Finiteness of the expression above follows from Propositions 7.1 and 7.3. Since both Propositions 
are very technical we postpone their statements to Section 7.1 and their proofs to Appendix D. □ 

Proof of Proposition 5.1. The Markov Renewal Theorem 5.5 yields 



lim — — 

t^oo ep[x) 



(xR > e s )ds = - 



g(u, v)X(dv)ir(du) 



S JR 



g{u, v)X(dv)ir(du) 



O0V 



ep(u) 



w' 3 ' 1 [P {uR >w) -NF {uCR > w)] dw — i—iridu) 



(uR > e") - NF (uCR > e v )] X(dv)n(du) 
1 



s Jo 



ep{u) 



1 



s Jo 



E 



w p - L [P (\uR\ > w) - NF (\uCR\ > w)] dw 



1 



N 



^^uCiRi + uQ 



N 



J2 [uQR 




ir(du) 



for all x G S. The penultimate line is justified by the symmetry of ep (Proposition 4.3), the last 
identity will again be justified by Propositions 7.1 and 7.3. 

With this identity, in view of [21, Lemma 9.3] we may easily unsmooth /, and one can finally 
infer that 

P 



lim e 0t l 

t—>00 



(xR > e 



2(31 



E 



N 



uCiRi + uQ 



i=l 



N 



i=l 



□ 



The remaining part of the proof of Theorem 2.7 and 2.9, i.e. the proof of the positivity of the 
limiting constant K is postponed to the Section 7. 

6. Similarities 

In the last two sections we considered general matrices and proved existence of the limits in 
Theorems 2.7 and 2.9. It still remains to prove nondegeneracy of those limits and the arguments 
will be given in Section 7. 

Now we consider similarities and our aim is to prove existence of the limit in Theorem 2.11. 
We assume now that the C, take their values in the similarities group and that the assumptions of 
Theorem 2.11 are satisfied. The idea of the proof resembles the previous case of general matrices, 
i.e. we reduce the problem to the renewal equation but this time we apply an extended version of 
the renewal theorem for random walks on G — M + x 0(d) (Theorem A.l [14]). 

We denote by ~jl the law of C. Since N J G ||</||' 3 7l(<i<7) = 1, the measure fip(dg) = N\\g\\P]Z(dg) is 
a probability measure and moreover rap — J G log UgH/i^df?) = A r E[||C|| ° log ||C||] > 0. By ~p* k we 
denote the fcth convolution power of 71, i.e. the law of Ilk- In the same way we introduce ^ k . Let 
Up denote the potential of /i p, i.e. Up = y,^l Hp ■ Then Up is a Radon measure on G and we 
have the following result: 



TAIL BEHAVIOUR OF SOLUTIONS OF FIXED POINT EQUATIONS 

Lemma 6.1. Let f g C c {M. d \ {0}), then for any a g G 

OO „ 

\\a\\^E[f(aR)] - (*„ * UpMf) = / Tpf(ag)4 k (d9), 

where 
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fe=o 



i>f(g) = lis 



- \\n\\-Pl 



N 



fi.gR) - E fiaCiRi 



i=l 



Proof. We write 
||o||-"E[/(ai2)] 



n-l 



N 



lair^E J2 (f(aL(v)R(v))-J2f(*L(vi)R(vi)) 
k=a L |„| =fc V j =1 

n-l p iV 

= ||o||-^JV fc E f(aU k R) - /(«H fc Qi?) 



E[f(aL(v)R)] 



J2 E[f(aL(v)R)] 



„ n-l r N i 

/ ^||ag|r% f{agR)-Y,f{*9CiR) bt N^dg) + \\a\\^ ]T E[/(aL(«)i?)] 
G fc=0 L i=l J k)|=n 



fc=0 



£ lj f (ag)rt k (dg) + \\a\\-e £ E[/(oL(«)iZ)] : 



Notice that for 5 such that m(s) < 1, assuming supp/ D B v (0) = 0, we have 



£ E[/(aL(«)iZ)] 
|v|=n 



< ciY«F[|an^| > „] < gffll < "WW 



(»?iioir t 



Cr? ]]«[[ -> 



Hence 



and so 



lim 

n— »oc 



£ E[/(oi(«)fl)] 

|?j|=n 



--8 - 



||o||-"E[/(aiJ)] =£ / M°9)l*?(dg) 



k=0 



□ 



To apply the renewal theorem we have to check that the function ijjf is directly Riemann inte- 
grable. We need the two following lemmas, whose proofs will be presented in Appendix C: 

Lemma 6.2. Let f g C^(M d \ {0}) and let n g N 6e fixed. Then there is a constant c = c(f) such 
that for every < e < 1 and every x\, . . ■ , x n , q £ K d 

/(xi + ...+:£„ + g) -^/(Xj) < c(^\q\ £ + ^ |xj| e |:rj| e j 



j'=i 



Lemma 6.3. Lei R be a solution of (1.2) smc/i i/ia£ /or every s < f3. E\R\ S < oo. Suppose that 
iVE||Ci|| /3 = 1, E|Q|' 3 < oo, then for any f g C^,(R d \ {0}), the function tpf is directly Riemann 
integrable on G = R + x K . 
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Proof of Theorem 2.11. By the renewal theorem (Theorem A.l, [14]), for any / G C c (R d \ {0}) 

hm ||o||-^E/(ail) = — / tl> f (g)dg, 

||a|K0;aSG t , TO/3 Jq 

where dg is the Haar measure on G normalized in such a way that for any radial function / on G: 
f G f(d)dg = f °°f(t)f. 

In the same way as in [14] we show that the convergence above is valid also for / <E Cc(K d \ {0}) 
and thus there exists a Radon measure A on W l \ {0} such that 

lim \\a\\- f) Ef(aR) = (f,A). 

\\a\\^0,a£G„ 

Since A is homogeneous it can be written in the form A = a <g> j^t, i.e. 

dt 



(/, A) := 



ii 



f(tw)a(dw) 



where a is a finite measure on S. 

Finally we have to justify the formula for a, i.e. to prove 

N 8 N 



a{S) = —E 



For this purpose take an arbitrary radial function / G Cc^R \ {0}). To simplify our notation 
define 

dt 



Then, on the one hand 



//«=/_ f(t)-r +s , 



dt 



(/, A) = / j fitwMdw)-^ = a(S) I f{t)j^ = a(S)I f (J3), 



dt 



and on the other, since ipf is also a radial function we have 

1 f°° dt 

(/,A} = - / 1>,®T- 
mp Jo t 

Thus 

(6.4) a(S)I f (fS) = — Tv/Wt- 

me Jo t 

Notice that for s < j3 we have 

J™ t-°Ef(t\R\) J =E\R\ S f(t)^ =m\'If(*)' 
Hence for s < /3, since E|_R| S < oo, we may write 



A 



ijji--X;iic,ito 



i=i 



AT 



2V 



i = l 



J/00- 



Now, letting s — > /? , in view of (6.4), we obtain 

1 Z" 00 r7/ 
(7(5)1/08) - lim — / t^ f (t)- 
«->-/3- mp Jo t 
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lim E 

s->/3- rap 



1 

] 



i=l i=l 
P N 



R, 



i A*) 



I AV- 



ID 



7. POSITIVITY OF K 

The aim of this section is to study positivity of the limiting constants. We have proven up to now 
that tails of solutions to equations (1.2) and (1.3) behave regularly at infinity, but we still do not 
know whether the limit is non-degenerate. Now we fill this gap. We start with some preliminary 
estimates, that will be used in the proofs. Next we consider similarities for which our results are 
much stronger and give a complete answer, i.e. we prove Propositions 2.12 and 2.13. Finally we 
consider general matrices and complete the proofs of Theorems 2.7 and 2.9. 

7.1. Moment bounds. Before we pass to proofs of positivity of the limiting constants we formulate 
here two Propositions containing useful estimates. However, since their proofs are long and technical 
we postpone them to Appendix D. 

Proposition 7.1. Assume that the following expectations are finite: 



supE 

yes 



(E0(s)) 
(El(s)) 
(E2(s)) 

Then for all y £ S , the identity 



N 



(7.2) 



2 



N 



yCiRi + yQ 



supE 

yes 



N 



Y^yCiRi + yQ 

N 

i=l 
N 

J2 \y°^ 



N 



supE 

yes 



i=l 
maxyCiR., 

i 

max yCiRi 

i 

t z ~ x (P (\yR\ >t)-NV (\yCR\ > t)) alt. 



holds true (for complex z = s+iv ). Both sides of the equation above are finite and define holomorphic 
functions in the infinite strip < Sftz < s provided that < s < s^ . Moreover 



sup 

yes Jo 



s-l 



>t)-NV (\yCR\ > t)\ dt < E0(s) + El(s) + E2(s). 



Proposition 7.3. For g < s^, there exists e > 0, such that E|i?| e 6 < oo implies finiteness of 
E0{g + e), El(g + s), E2(g + e). 

7.2. Similarities. 

Proof of Proposition 2.12. First note, that K + > if and only if u{S) > 0. We are going to use 
Propositions 7.1 and 7.3 in the formulation adjusted to the similarities. Namely, let 



(E0(s)') 



supE 

yes 



N 



C > R > + Q 



N 



Y2 CiRi 



1=1 
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N 



(El(s)') 



(E2(s)') 



supE 

yes 



supE 

yes 



E CiRi 

i=l 
N 



— max \ CjRj 



max \ CiRi 



Then Proposition 7.3 holds with the same proof, and the finiteness of (E0(s)'), (El(s)'), (E2(s)') 
implies, as before, the identity 



(7.4) 



E 



'E CiRi + q 



N \ c 

Eiw ^ / 

i=l / J o 



i 2 " 1 (P (|i?| >t)-NP {\CR\ > t)) dt, 



for complex z = s + iv\ moreover, both sides are holomorphic for < f3 + e. On the other hand, 
we have for complex z = s + iv, a < s < (3, since E |i?| s < oo, 



E 



N 



\R\ 



E 



= E 



A' 



\R\'-^\Ci 



\Ri 



= (1 -m{z))E\R\ 2 



m\ z -Eli\CiR 



E\R\ Z 



for 5Rz < /?. 



1 — m(z) 

Now the rest of the proof of Proposition 2.12 is the same as the arguments in [14], but we include 
it here for completeness. 

Suppose that K = 0. Then both the numerator and the denominator of the right hand side are 
holomorphic in < SRz < ft + e and 1 — m(z) has a simple zero at (3 (m'(/3) ^ 0). Therefore, the 
right hand side is holomorphic for 5Rz < /3 + e. On the other hand, if z — s € R, the left hand side 
is the Mellin transform j(s) of the law of \R\, which is well defined for s < 6^ called the abscissa 
of convergence of j. The Landau theorem (see [36], Theorem 5, page 57) says that 7 cannot be 
extended holomorphically to a neighborhood of 9oo. Hence ft + e < #oo, i.e. E li?^ 6 < 00. Now 
suppose that 6*00 < s^. Then (7.4) is holomorphic for Sftz < 6*00 + e, and so repeating the above 



argument (using that m(z) 7^ for z ^ a, (3) we get a contradiction. 

7.3. A sufficient condition for the limit to be positive. We introduce the r.v. 

jv 

B:=J2 CiRi + Q> 



□ 



since we will use some features of the afhne equation R = C\R + B in the sequel. 

Lemma 7.5. // (not triviality) holds, then P (dr + B = r) < 1, for all r e R d . 

Proof. Condition (not triviality) states that no Dirac measure solves the fixed point equation, so R 
takes at least two different values. If now 

N 



C ir -Q = J2 G i R i 



one sees by conditioning on (Ci, . . . , Cat, Q) (then the LHS is constant), that this identity can 
only hold on a set where all values of the Ri are fixed. But, since they are independent and non 
degenerated, this set has probability smaller than one. So P (C\r + B = r) < 1 for all r € M. d . □ 
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Proof of Proposition 2.13. Obviously, if (not triviality) does not hold, the solution is a constant, and 
a(S) = 0. Conversely, assume that a(S) = 0, so by Proposition 2.12, E|_R| i < oo for all s < Soo- 
With the definition of B, R satisfies the affine stochastic fixed point equation 

R = CiR + B, 

which was studied in [14]. By the assumption E|Q| 7 < oo, we also have E|5| 7 < oo. But then 
(not triviality) may not hold, since otherwise condition H of [14] is satisfied, and it is shown in [12, 
Proposition 2.6], that then E |i?| 7 = oo, which would be a contradiction. □ 



7.4. General matrices. In the case of general matrices, the crucial identity becomes more subtle. 
Using T s v s = ^jp-P s and E |i?| s < oo, m(s) < 1 for all s e (a,min{/3, Soo}) we have 

m{s) [ E \yR\ s v s {dy) = N [ T s E\yR\ s v s (dy) =V / E \yC t Ri\ s v s (dy), 
Js Js i=1 Js 

and thus the identity 

(7.6) jf E \yR\ s v s {dy) = (1 - m(s))~ 1 J^E (^\yR\ s - ^ \yCiIuA v s {dy). 

We still can show that the right hand side has a holomorphic extension around j3 if and only if K = 
0, but unfortunately, we cannot use the Landau lemma, because the LHS is not a Mellin transform, 
only a mixture of those. Nevertheless, what we can show by this argument is: If E = oo, then 
the right hand side cannot have a holomorphic extension, thus K > 0. (Note that E \R\ < oo 
readily implies K = 0.) 

Proposition 7.7. There is S > 0, such that 

N 



K(z) := Je \^yR\ z -Y^yCiRiA u z (dy), 



is a holomorphic function on Bs{/3). 

Proof. By Proposition 7.3, there is e > such that E0((3 + e),El{/3 + e), E2(fi + e) are finite. By 
Proposition 7.1 then for each y g 5, 



k{y,z) :=E 



v 



^ yC t R t + yQ 



N 



J2 \y°* R * 



is holomorphic. By Corollary 4.12, there is a (possibly smaller) S > 0, such that v z is a holomorphic 
function on B$((3). But then 



K(z) = J k(y,z)v z {dy), 

also defines a holomorphic function on Bs(f3). □ 

Also by Corollary 4.12 , m(s) = Nk(s) is holomorphic on B$((3), with m(/3) = 1. Since m'(/3) > 
(by convexity), the right hand side in (7.6) has a holomorphic extension around j3 if and only if K = 0. 
Assume that K = 0. The left hand side is a holomorphic function on Bg((3) n {5Rz < /?}, which is a 
domain in C, thus a holomorphic extension of the right hand side is also a holomorphic extension of 
the left hand side. Let us call it Unfortunately we cannot conclude that ^f(s) = J S E \yR\ s v s (dy) 
for s > P, this is caused by the extra dependence on z appearing in v z and this is why the Landau 
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lemma cannot be used here. We only know, that \&(s) = J s E \yR\ s v s (dy) for s < (3. But that at 
least implies that if ^ exists on Bs(j3), then 



lim / E \yR\ s v s (dy) = lim \Ks) < oo. 

«t/3 Js s^/ 3 



So the following lemma gives the contradiction if E \R\ = oo, implying that K > at the outset. 
Lemma 7.8. IfE\Rf = oo, then 

(7.9) lim / E\yR\ s v s (dy) = oo. 

Js 

Proof. Now if E i?)"" = oo, then there is xo £ 5 with E |xoi?| CT = oo. By Lemma 7.10, then there is 
already e > such that E \yR\ a = oo for all y £ B £ (xo) O S. By Lemma 4.4, 

liminf / E \yR\ s v s (dy) > liminf / E (|?/i?| s l/i „fl|>u ) v 8 {dy) 
s^a J s s^a J s llw >- 

> liminf p / E (\yR\ s l{k/K|>i}) <t>{dy) 
s -> a Js 

>p / liminf E (\yR\ s 1{\ v r\>i}) 4>(dy) = oo. 



In the penultimate line, we used Fatou's lemma. □ 

Lemma 7.10. 7/E (\xR\ s ) — oo for some x S S, then there exists e > such that E (lyi?^) = oo 
for all y G B E {x) . 

Proof. Case 1: Suppose that there is an open U C R d such that E |yi?| s = oo for all y £ U \ {0}. 



Then of course E 

Case 2: If not, then 



\y\ 



= oo and the image of U \ {0} on the sphere S is open. 



B := {y eR d : E\yR\ s < oo}, 

is dense. But then there is already a basis yi, . . . , yd of M d such that E \yiR\ s < oo for all 1 < i < d, 
and so E \yR\ s < oo for every y £ R d , but this is contradiction. To construct such a basis just note 
that the complement of the linear hull lin(j/i, . . . , yk) is a hyperplane, so we can always find a vector 
in B which is independent of the previously chosen basis vectors. 

□ 



8. Applications of the generalized smoothing transform 
8.1. Multivariate examples from kinetic gas theory. 

Particle velocity in a homogeneous gas. Consider the distribution V £ M. d of particle velocities in a 
homogeneous gas. If two particles collide, the interacting particles change their speeds vi and v- 2 to 
the post-collision speeds v[ and v' 2 according to the formula 

v[ = C^vi + C^v 2 , v' 2 = C[ 2) v 2 + C^V, 

where (C[ 1 \C 2 1 ' > ) and (C[ 2 ^ , C 2 2 ^ ) are d x d-matrices, transmitting the interaction and depending 
on the scattering angle. To obtain a stochastic fixed point equation, we interpret the interacting 
particles as randomly picked, so we may assume that v\ and v 2 are iid and distributed according 
to V, and independent of the scattering angle, say the matrices, which can also be assumed to 
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be random. So if (Cijd) denotes a generic copy of (C^^C^), in the equilibrium the speed 
distribution V would satisfy the stochastic fixed point equation 

V = C1V1 +C 2 V 2 , 
with V, V\, V 2 iid and independent of (C\, C 2 ). 

A particular model: Maxwell-type gas. Bassetti and Matth.es [9, Example 6.1] study the following 
model for a Maxwell-type gas: 

Ci = UY T Y, C 2 = Id - UY T Y, 

where Y is a random unit row vector, modelling the line of collision and U G Rj!~ is a random 
variable modelling the inelasticity. They assume <f> = C(Y) to be the uniform distribution on 5*, and 
independent of U. Moreover, they assume that 

E[U(l -U)]=0 and E[U 2 (l - U) 2 } < E[\(w u Y)\ 4 ]EU 2 , 

where w\ denotes the unit vector (1, 0, . . . , 0). Let's finally assume that the distribution ip of log U 
is nonarithmetic. Then we have that for all x G S, 

P((iC0~ e -.loglxCil e •) = <£«>V>, 

and since C{C) = |£(d) + |£(C 2 ), 

V KeS P ((sCT e ■, log G •) > ^ ® V, 

thus (minorization) holds. Since is the uniform distribution Z on S 1 , also (not vanishing) holds. The 
same would be true for any distribution (f> having a density w.r.t. to I. 

By the very definition of Ci, C 2 , S = Id satisfies (covariance). Now lets consider m. We have 

m(2) =2E[sup |Cx| 2 ] = E[sup((C x a;, da) + (C 2 i,C 2 i))] = E[supa; T (C 1 T Ci + CjC 2 )x] 

x£S x&S x£S 

<sup.x T Idx-2E[[7(l - t/)]E[inf xY T Yx T ] = x!dx T = (x,x) = 1. 

Moreover, it follows from [9, 6.1.2], that m(4) < 1. Thus a = 2, and the assumptions of Proposition 
2.4 are satisfied, giving a unique solution V. 

In [9, 6.1.4], it is shown that m(s) > 1 for some s > 4. This gives the existence of j3 with 
m(/3) = l,m'(/3) > for some /? > 4. Thus our Theorem 2.9 is applicable in this situation. 

Precise tail estimates for the Maxwell-type gas. 

Lemma 8.1. In the situation above, Soo = sup{s > : Et/ S < 00}, and if I denotes the uniform 
distribution on the sphere, 

e s = 1, v s = I V.s < Soo, 

as well as 

m(s) = E (U s + \wi - UY\ S ) . 

Proof. We have 

E ||C|r = ^E Hd II s + ^E ||C 2 |r < \w s + 1(1 + EU S ) = \+ W\ 

and also E \\C\\ S > ±E ||Ci|| s = ±EU S . 

Since Y has the uniform distribution I on S, we infer that for all x G S 

C(x - xUY T Y) = C(wi - UY), C(xUY T Y) = C{UY), 
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thus for all s < Soo, / G C(S), 

T s f(x) =E (f((xCD \xC\ s ) = iffi (fdxC.D \xC x \ s ) + ±E (/((a:C a )~) \xC 2 \ s ) 

=±E (/((F)-) If/D + l -E (f(( Wl UYT) \ Wl - UY\ S ) . 

Since the RHS does not depend on x, e s = and v s = I for all s < Soo. 

The formula for m(s) follows from Proposition 4.3, using that T s lg = k{s) and m(s) = 2k(s). □ 

Given the distribution of U, /3 is explicitly computable from this formula. 
Finally, in this particular example, identity 7.6 reduces to 

/ E \yV\ s u s (dy) = E |Yy| s = (1 - mfc))"^ (|Y(CiVi + C 2 V 2 )\ S - |YCiVi|' - \YC 2 V 2 \ S ) , 

with V\,V 2 iid and distributed like the solution V. Thus the LHS is the Mellin transform of the 
random variable |YV|, and our holomorphic argument works as in the situation of similarities. Since 
by [9, 6.1.4], there exists s < Sqo such that E \V\ S = oo, we deduce that necessarily E IV] 13 = oo. 
Thus an application of Theorem 2.9 yields that 

lim i^P (xV > t) = K > VieS. 

8.1.1. A second example with similarities. In the same paper by Bassetti and Matthes, there is also 
an example that considers the multivariate smoothing transform with similarities, [9, Example 6.2]: 
It is stated as follows: Let A, B random matrices in SO(M. d ), and a,b be non-negative random 
variables, independent of (A,B) and satisfying 

m(2) =E[a 2 + b 2 ] = 1, m{s) = E[a s + b s ] < 1 

for some 2 < s < 3. Then they consider the homogeneous smoothing transform associated with 

Ci = aA, C 2 = bB. 

One readily checks that, since A, B are orthogonal, 

E[dldCj + C 2 IdC 2 T ] = E[a 2 AA T + b' 2 BB T ] = E[a 2 + b 2 Id] = Id. 

Thus by Proposition 2.4 there is a unique solution V, and if there is j3 > 2 such that E[a^ + b /3 ] = 1, 
then by Theorem 2.11, V has heavy tails with tail index j3. 

8.2. The inhomogeneous smoothing transform with (3 = 2. Consider the inhomogeneous 
fixed point equation 

N 

(SFPE) R^Y^CuRk + Q. 

fe=i 

Let EQ = 0, < Var Q < oo, and assume f3 = 2. Then by Proposition 2.5, there is a unique solution 
R with ER = 0. Notice that this R = in the case of the associated homogeneous equation. Also, 
R = Q = would be the only possible solution when only considering solutions on the positive cone 
as in [13, 30]. But in our model, R is certainly nondegenerate due to the influence of Q. 

This is a very nice example to study, because another feature, which is only possible in the 
multidimensional setting, occurs: If one would allow for noninvertible matrices C, and 

[supp C] ■ (supp Q) = {0}, 

where [supp C] denotes the smallest closed semigroup generated by supp C, and (supp Q) the linear 

hull of supp Q; then the solution to (1.2) is given by R = Q. 
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In our model, we only consider regular matrices, so this situation will not occur. Nevertheless, 
one may expect that properties of R strongly depend on the interaction between the distributions 
of C and Q. In the case = 2, this can be made explicit (by the way showing, that K is positive 
in this case): 

Lemma 8.2. Consider the inhomogeneous SFPE (1.2). Let the assumptions of Theorem 2.7 hold, 
and assume EQ = 0, VarQ > and (3 = 2. Then the unique solution R with MR = satisfies 

e(x) 



lim t 2 P (xR > t) 



4/, 



E(yQy V2 (dy) > 0, 



for all x £ S . 



Proof. We calculate for y G S, using the independence of Ri, . . . , Rn and (Ci, . . . , Cat, Q) as well 
as ER = : 

2 N 

2 



N 



e J2 y c i R i + vQ 



- J2 \yCiRi 



i=i 

N 



=e ^(yCAXyQi?,) + J2(y c * R *)vQ + (yQ)' 



N 



= E (E [ {yCiRi ) (yC, i^) | Q, Cj}) + E (E [ (yC^ ) (yQ) | Ci , Q] ) + E{yQf 



i=l 



N 



= ^2E((yC i -ER)(yC j m)) + (( y C t ER)(yQ)) + E(yQ) 2 
=E(yQ) 2 . 

Now by Proposition 5.1, the limit is given by J s E(yQ) 2 ^(dy) . It is shown in Lemma 4.4, that 
supp(^) = S, thus the last integral is indeed positive. □ 

To interpret this formula, remember that v 2 is the invariant measure of the operator T2 defined 

by 

T 2 f{x)=E(\xC\ 2 f{{xCT)). 
So it is directly connected to the action of C. 

8.3. On contraction conditions. In their survey article [33], Neininger and Ruschendorf intro- 
duce several probability metrics, namely Zolotarev and minimal Z p -metric, and discuss the resulting 
conditions for S to be a contraction on *Xt2 (0). 

They list several conditions for 5* to have a unique fixed point with zero expectation and finite 
variance, namely 

(!) E Efc=i llCfcll < !» which is the 

one we used and comes from the Zolotarev metric, 
(2) X^fcli E ll^/T^H < 1' wm ch comes from the minimal ^-metric, and 



(3) Ef =1 ECjC 



< 1. 



The last one is the weakest. In the context of the contraction method (see there for details), they 
pose the question, which is the best, i.e. weakest condition possible. See [33, Problem 3.2]. 

We cannot resolve their question, but our example above gives a good hint, which may be the 
optimal condition: We suppose it is m(2) < 1; for the example above shows, that if m(2) = 1, the 
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unique solution with zero expectation has infinite variance. In particular, if m(2) = 1, there is no 
solution with finite variance. 



Appendix A. Reduction to identically distributed weights 



Proposition A.l. Let a be a r.v. with uniform distribution on the symmetric group of order N , 
independent of all other occuring r.v.s. Then R (with iid copies Ri, . . . ,Rn) is a solution to (1.2) 
if and only if 



N 



(A.2) 



; = 1 



i.e. R solves the SFPE associated with the vector (C CT (i), . . . , C a ^), Q)- 

Moreover, r £ M. d satisfies (eigenvalue) for (Ci, . . . , Cjv, Q) if and only if r satisfies (eigenvalue) 
for (C ail)l ...,C a(N)l Q). 

Proof. Let R\, . . . , Rn be iid, independent of (Ci, . . . ,Cn,Q) and a as defined above. Then for all 



N 



Ef J2 °k R k + Q 



E 



E 



= E 



Ef (j2 C* {k)Ra{k) + Q 

E^f(^pC a(k) R k + Q 



(Ci, . . . , Cn, Q), (T 
(Ci, . . . ,C N ,Q),a 



Now the distribution of R is a solution to the original equation, if Ef(R) equals the first line, and 
it is a solution to the permuted equation, if Ef(R) equals the last line; so both equations are indeed 
equivalent. 

The statement about the eigenvalues can be readily checked by taking expectations in both 
equations. □ 

Remark A.l. Proposition A.l shows that we may w.l.o.g. assume that C\, . . . , Cn are identically 
distributed (but dependent) with generic copy C. 



Appendix B. Existence and Uniqueness of Solutions 

The aim of this section is to show existence and uniqueness of solutions to (1.2) and (1.3), given 
by Propositions 2.4 and 2.5. For this purpose we consider the mapping S as a contraction operator 
on an appropriate complete metric space consisting of probability measures. Then both existence 
and uniqueness of solutions follow from the Banach fixed point theorem. 

B.l. The Zolotarev metric. The right metric on measures for our purpose is the Zolotarev metric 
Cs, introduced by Zolotarev [37], and considered in this context by Rachev and Ruschendorff [34]. 
For X, Y two random variables on a common probability space, we define 



(B.l) 



(s(Z,Y) '■= sup{|E (f(X) - f(Y))\ : 
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where D s is the space of functions defined as follows 

(B.2) V s = {feC k (R d ) : V^ eRd \f w (x)-f {k Hv)\<C\x-y\'- k } with fc = Ts - l] , 

and \a] denotes the smallest integer > a. Observe that for < s < 1, D s is the set of s- Holder 
functions on M. d . 

For properties of the Zolotarev metric £ s see [32, Section 2]; or for a recent treatment of the 
Zolotarev metric on general separable Hilbert spaces see [19]. It turns out that on an appropriately 
defined measure space the Zolotarev metric is complete and the mapping S is a contraction. The 
details are as follows. 

Let a = (oi, . . . , ad) 6 N|j be a d-dimensional multi-index, and set |a| = ai + • • • + a^. For a 
probability measure v on R d , denote by M k (v) the sequence of mixed integer moments up to order 
k, i.e. 

M k {v) = 1 / x\ x ■■■x a d d dv(xi,...,x d ) I 

\J«- d / |a|<fc 



For a fixed mixed moment sequence M k , we introduce the following subsets of the space of 
probability measures on M. d (as before, k = \s — 1]): 

(B.3) m d s {M k ) := L : f \x\ s v(dx) < oo, M k (v) = M k 



Thus, we consider subspaces of measures such that we fix integer moments, i.e. the expectation if 
1 < s < 2, and also the covariance matrix if 2 < s < 3, and so on. 
Then we have the following Lemma: 

Lemma B.4. The metric spaces (dJl d (M k ), £,) are complete for s > 0. Moreover, ifm(s) < 1, then 
for some uq the operator S n ° is a contraction on the space (9Jtf {Mk), Cs)- 

Proof. The first part of the Lemma was indeed proved in [19, Theorem 5.1]. However for reader's 
convenience we give some ideas of the proof. 

Let v n be a Cauclry sequence in 9Jl d (M k ) with respect to the Zolotarev metric. By [37, Theorem 5], 
it is also a Cauchy sequence in the Prokhorov metric, thus weakly convergent towards a probability 
measure v - observe that characteristic functions belong to J) s for any s > 0. By [37, Theorem 6], the 
sequence of J \\x\\ s v n {dx) is also a (real) Cauchy sequence, thus bounded, implying / [|a;|| s v(dx) < 
oo. For s > 1, this also yields uniform integrability of the moments of order up to s, thus v again is 
in DJl d (M k ). It remains to prove that Cs{v n , v) — > 0, for which we refer to [19, Lemma 5.6]. 

The second part follows directly from the inequality proved in [32, Lemma 3.1]: 
C(Sv,Sr,) < ^E^||C i || s jC s (i/^) = iVE||C|rC s (^r?)- 

□ 

B.2. Proof of Existence and Uniqueness. Now we are ready to prove both Propositions. 

Proof of Propositions 2.4 and 2.5. If m(s) < 1 for s < 2, the mapping S n ° is a contraction on SDtf, 
9Jtj(r) resp. 9Jl d (r, S) which is a complete metric space (Lemma B.4). Therefore by the Banach fixed 
point theorem there exists a unique fixed point R of 5, being a unique solution of (1.2), resp. (1.3). 
Notice that for 1 < s < 2 (2 < s < 3), we have to fix the expectation (and the covariance matrix) in 
order for iS to be a contraction. Moreover, to ensure that S maps 9Jl~(r) into 9Jtj(r) (9Jtf (r, E) into 
9Jtj(r,S)), we have to know that r satisfies condition (eigenvalue) (S satisfies (covariance) for the 
homogeneous equation). Observe further that in the homogeneous case R = is a trivial solution 
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of (1.3). Since for s < 1, So is an element of 9Jtf, there are no other solutions of this equation. Thus 
(1.3) possesses nontrivial solutions with finite expectation only if a > 1 and thus s > 1. 

To prove the second part of both Propositions, observe we already know that all moments of R 
up to some sq > a are finite. By convexity of to, to(s) < 1 for all s g (sq, min{/3, Soo}). Then 
arguing as above, Lemma B.4 assures the existence of a solution with finite s-moments, which by 
uniqueness coincides with R. 

If s > 2, then S is a contraction only on subspaces of measures with fixed mixed integer moments 
of higher order; and we have to check whether S is still a self-mapping of these subspaces. But this 
is no problem, since Lemma B.4 applied to p = \_s\, (here and subsequently [_a\ denotes the biggest 
integer < a) gives finiteness of all mixed moments of R up to order p, and since R is a solution, 
then also S(M p (R)) = M p {R). Starting with p = 2, we can proceed in the same manner at every 



integer step. 

Appendix C. Proofs of Lemmas 6.2 and 6.3 
Proof of Lemma 6.2. We estimate separately 



□ 



(C.l) 
and 



f{xx + . . . +x n + q) - f(xi + . . . + x n )\ < c\q\ e 
f(xi + . . 



(C.2) /(si + .••+*») <cX>i| £ |z;| £ - 

3=1 i^j 

On one hand (C.l) is bounded and on the other it is dominated by c\q\. Hence 

+...+x n + q)- f( Xl + ... + x n )\ < c(f)\q\ e . 

For (C.2), we first prove 

(C3) \.f(x 1 +x 2 )-f(x 1 )-f(x 2 )\<c(f)\x 1 \ £ \x 2 \" 
and then we proceed by induction as follows 



f(x 1 + . . . + x n ) - }] /(a 

3 = 1 



< \ f(xi + ... + !„)- f(xi) - f(x 2 + ■■■ + x n ) \ + 

< c|a;i | e |a; 2 + • • • + x n \ e 4 



E 



-£/(* 

3=2 



i#3';i,3>2 



For (C.3) we write w = (ui, . . . , u^), y = : (j/i, ■ • ■ , 2/d) £ K d and 



f(u + y)-f(u) 

m 



(Is 

M - /(o) 



/(m + sy)ds 



»1 d 

/ y2 d jf( u + s y)yj ds i 

Jo , = i 



1 d 



^2djf{sy)yjds. 



3 = 1 



Now notice that 



d 3 f{u + sy) -djf(sy) 



1 d 



^ d k djf(ru + sy)yjU k dr. 



k=l 
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Hence 



\f(u + y)-f(u)-f(y)\< [ [ V Idkdjfiru + sy^UkWyjldsdrKCluWy] 

On the other hand \f(u + y) — f(u) — f(y)\ is bounded, so (C.3) follows. 
In several places, we will use the inequality 



□ 



(C.4) 



valid for yi > and s > 0. For s < 1, this follows from the subadditivity of x h-> x s . For s > 1, 
it follows from Jensen inequality, when applied to the discrete probability measure — Y^i=i ^vt an( ^ 
the convex function x t— > x s . 

Proof of Lemma 6.3. Writing g = tk, t € K. + , k <G K we have to show that 

oo 

£ sup \tpf(tk)\ < oo. 

„ = _ 00 e»<t<e»+»;fceJC 

Assuming supp/ C K d \ -Br/(0) and e™ < t < e n+1 by Lemma 6.2 we have 



|V>/(ifc)| < t- p I 
< ct~ p E 



N 



v 



£ ItkCiRifltkCjRtf + \tkQ\ E ) • IfEICi^i+iQi^t-i} 
< ce -(^)«£E[|Qi^|C^f l^lC^I+IQI^e— i}] 

+ce-^- £ )"E[|Q| £ l {y: | Ctfl( | + |Q| >lje -»-i } ] . 

Denote the first expression above by I(n) and the second one by II{n). Let the random variable no 
be defined as no = [log?/ — 1 — log(^ ■ |C;i?i| + |Q|)] . The value of a constant c. may change from 
line to line, but is finite and depends only on the indicated variables. Then 



< cE 



\Q\ S - £ 



-(/3-e)r 



n>no 



< cE 



\Q\ e e 



e-(/3-e)n 



< On 



V ,0-e 

J2\CiRi\ + \Q\) \Q\ £ 

• i=l ' 

^ i=l ' 

where the last line follows by an application of (C.4). The last expression is finite since 

E[||C i ||' 3 - £ | J R/- £ |Q| e ] =E[| J R i |^- £ ]E[||C i ||^- e |Q| £ ] < E[|i?/- £ ] (E\\df ) V (e|Q|") ; 
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Next, recalling that i ^ j and Ci are identically distributed we write 



< Cat! 



Y, e-^- 2 ^ n \C x R x \ e \C 2 R 2 \ 



n>nn 



< 



c A rE[e-(' 3 - 2e )" |C 1 i? 1 HC , 2 i? 2 



< c N „E 



N 



k=l 



P-2e 



< c Nn N w ^~ 2 ^ 



< 



< 



J2\ C k R k\ + \Q\) \CiRi\ e \C 2 R 2 \ £ 

N . 

]T \C k R k f- 2E + \Q\ f3 -' 2E ) |C 1 i? 1 nc 2 i? 2 | 

\Qf- ae \C 1 R 1 \ B \C 2 R2\ 



k=l 



r N 



Y\CkRk\ - 2e \CxRxnC 2 R 2 
■fc=i 

Cj v i ^[(e| j R| £ ) 2 (e|Q|^^(e||C 1 ||^ 
+ ^(E||C fc f) " (E||Cif) 7r .E[|i? fc |''- 26 |iii| s |iJ 



fc=i 

where we used again (C.4), and the general Holder inequality in the last line. Only finiteness of 
the last term requires some proof. If k ^ {1,2}, then 

E^Ruf-^R^Rtf] = EIJ?* I^-^EI^! | E E| J2 2 | e , 

otherwise 

EftRkf-^lRiFW] = E\R^- s m\R 2 \ e , 
and both expressions are finite. □ 

Appendix D. Proofs of Propositions 7.1 and 7.3 
To prove Proposition 7.1 we need the following lemma: 
Lemma D.l. For all < s < Soo, 



r 



(max\yCiRi\ > tj -N¥{\yCR\ > t) dt < (E2(s)). 



If (E2(s)) is finite, then 
(D.2) 



z j t z ~ x (P (max \yCiRi\ > tj - N¥ {\yCR\ > tfj dt = E 



N 



max yCiRi 



\y°^ 



i=l 



for all z S C such that < < s. 

Proof. If (E2(s)) is finite, then the function 



N 



"-ME 1 



{\yC,R,>t\} — - L {max i |yCifl,|>t} j 
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2') 



is also absolutely integrable with respect to A ( 



E 



N 



't* 1 ( X] ^{\yCiR,\>t} - l{maxi|»CjRi|>t} 
\i=l J 



(It 



1*1 



{|2/C ! fl ! |>t} - l{max i |j/C i fl i |>t} ] dt = \z\ • (E2(s)), 



Thus Fubini theorem applies and we obtain 

z t z ~ l (p (max\yCiRi\ > tj -N¥(\yCR\ > t)j dt 

= jf zt z ~ x (f (max \yCiRi\ > tj - ^ P {\yCiRi\ > i) j dt 



( N \ 

zt Z L E l{max l \yC,R,\>t} ~ ^ 1 {\yC l R i \>t} I * 



=E 



maxyC,i?, 



A' 



i=l 
n /.Ij/Ciflil 



zt z ~ l dt 



Proof of Proposition 7. 1 . We have 



\zt z 



N 



<s / t 
10 



s-1 



S I t 

'o 

5 / 

10 



s-1 



i=l 
N 



i=l 



- JVP(|vCJ2| >i) 
v 



dt 



dt 



> t 



r/f 



s-1 



(max lyC.i?,! > tj - NV {\yCR\ > t) 



dt 



□ 



Now, in view of [21, Lemma 9.4] the first expression is bounded by (E0(s)) and by Lemma D.l the 
second and the third one by (El(s)) and (E2(s)), respectively. 

The same calculation without absolute value signs yields the identity. Turning to the last asser- 
tion, since = |z|t s_1 , for all z with ¥tz < s 



zt z 



(\yR\ > t) - ATP (\yCR\ >t)\dt<D, 



for some finite constant D > 0. 
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Now if 7 is some closed curve in the domain < < s, then we may change the order of 
integrals and obtain: 

r r°° 

zt 7 - 1 (P (\yR\ >t)-NV (\yCR\ > t)) dtdz 



7 JO 

J zt z ~ l dz^J (P (\yR\ > t) — NP (\yCR\ > t)) dt = 0, 

and holomorphicity follows. □ 

Proposition 7.3 follow from the following lemmas 
Lemma D.3. For s < Soo , E0(s) is finite if: E |i?| OV(s_1) < oo. 
Lemma D.4. For s < s^, E2(s) is finite if: There is s/2 < 7 < s s.t. E |i?| 7 < 00. 
Lemma D.5. For s < Sqo, E1(s) is finite if: 
Case a < 2: E\R\^ < 00. 

Case s > 2: E |i?| 7 < 00 for 7 := max{s — 1,8- j^2\} = s — 1. 

Unfortunately, proofs of all lemmas are long, therefore we first present how they imply the 
Proposition and later we give their proofs. 

Proof of Proposition 7.3. The choice of e is described by the previous lemmata. Of course g+e < s,^, 
the other restrictions can be reduced to: 



(1) g — e > g + e — 1 (Lemma D.3, Lemma D.5 for g + e > 2) 

(2) g — e > i(g + e) (Lemma D.4, Lemma D.5 for g + e < 2) 

These restrictions can always be fulfilled. 

Proof of Lemma D.3. If s < 1, by the inequality \a s — b s \ < \a — b\ s , valid for a, b > 0, 



□ 



E 



N 



yCiRi + yQ 



1=1 



< E\Q\ S < 00. 



If s > 1, by the inequality \a s — b s \ < s \ a — b\ max{a s , b s }, valid for a, b > 0, 



E 



N 



^yC.R. + yQ 



i=l 



iV 



^yQRi 



i=l 



< sE |Q|max 



N 



yC t R t + yQ 



s-l 



N 



i=l 



N 



< sE mqi ^2 + \Q\J 

< S (n + i) ov ( s - 2 >e (^2 iCiikr 1 101 + ion 

< s^ + l^^Eli?! 5 '-^ f^HCill'- 1 ^ +E|Q| 4 

< s(Ar + l) 0v ( s - 2 )E|i?| s " 1 iV(E||C l |r) £ ^(E|Qni + E| 



In *, we used subadditivity of x n- x s 1 for s < 2, and (C.4) for s > 2. 



□ 
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Proof of Lemma D.J^. Let J- := <r(Ci, . . . , Cn). Observe that the function is nonnegative, so we 
may change integrals: 



v 



> t) ] f-^-dt 



(E2(s)) =E \J2 lw<W - ^max^ \yC.R, 

=s J Q (j2 v (iw^i > *) - p (jggSt \y CiRi 

=sE ^jf° f^P QyCiRi] > t\F) - P ( m« \yGAl > t\F)J t^dt^j 

(D.6) =sE ^ ^P ^maxJyC-i^l < t\fj - 1 + ^ P {\yC l R l \ > t\ J") J ^d/j 

Conditioned on J 7 , the r.v.s |yCji?j| are independent, so 

, x v 

max li/CiiZil < t\ T = TT (1 - P (lyC.i?,! > i\T)) 

l<i<N j A A 



Exactly as in [26, proof of Lemma 4.7], we use the inequality 1 — x < e 21 , valid for x > 0, to obtain 



N 



? =1 V(\yCiRi\>t\T) 



7 



i=i 

For s/2 < 7 < s, such that E |i?| 7 < oo, the Markov inequality yields 

N N N 

J2^(\yCiRi\ >t\F) <Y,^{\yCiRi?\^)^ <t-m\R\^J2\\ c i 

i—1 i—1 i—1 

Define the function g(x) := e~ x — 1 + x and note that g(x) is increasing for x > 0. Thus 

N N 

9 (J2^(\yCiRi\ > t\?)) < g(rm\RrJ2\\c*\r)- 

i=l i=l 

With the change of variables u = t">E |i?| 7 J2i=i W^W* , the 

inner integral in (D.6) can be estimated 

by 

(E |i?| 7 NE \\Cif) sh / (e- u - 1 + u^-^^du. 



The integral is finite ifl < s/7 < 2, which was exactly our choice. Finally, evaluating the 
expectation in (D.6), we see that for some finite constant c > 



(E2(s)) < cN s ^ (ElRpy^EiWdf)^ < c (E |i?| 7 ) 5 ' /7 N S ^E (\\C\\ S ) . 



□ 



Proof of Lemma D.5. We will only consider H(s), defined below since (El(s)) < H(s) + (E2(s)). 



H(s) :=E 



N 



N 
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N 



Y \yCiRi\ 2 + Y yCkRkyQRi - E \yCiRi 



3 



k=tl 

i n 



f N \ 2 N 

Ely^l 2 -Y^^yCiRi 



2\ 2 



\i=l 



=: H 1 +H 2 . 



Case s > 2: We have 



N 



Hi | E \yCkR k \ \yCiRi\ Y \y c A\ 2 + Y \yCkR k yQRi 

k=jtl \i=l k^l 



E 



t=l 

-1 



< s2 *- 2 iv4- 2 E [ (e i*w a ) f E iw^^i 

\\i=l / \k=£l 



s2i~ 2 (N 2 - N)i-'E ( E \yC k R k \ § \yQRt\ 
k^i 



In the first line, we used the mean value theorem, in the second and third line (C.4). Now by 
application of the (generalised) Holder inequality, this expression is bounded by terms of E |i?| s , 
E|i?|,E||C|| s . 

To estimate Hi we adapt arguments from [26, Lemma 4.1]). Set q = s/2, p = \q] and 7 = q/p £ 
(j)/p+ 1, 1]. Define the set of multiindices 

A p (N):={( Jll ...J N )eN N : j x + ■ ■ ■ + j N = p; < j; < p - 1}. 



For any sequence j/,- > 0, 
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P 
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y 3 i ■■■ y 3 N 
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Using this inequality and the conditional Jensen inequality, we obtain 
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e Eiy^n ~Y\y CiR i 
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Case s < 2: For by the usual inequality \a s / 2 — W 2 | < \a — 6| s ^ 2 and the Cauchy-Schwartz 
inequality: 

m KEY^lyCkRk^ lyCtR^ < N 2 (E\\C\\ S ) (e\r\ 

k^l 

Finally 



H 2 <e iY^\yCiRi\ s 

\i=l 

( N 

<2E lydXi 



max \yCiRj 
i<i<N ly 



N 



m^lvdRil ) ) <2(E2(s)). 



□ 
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